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We explore the effects of Coulomb interaction upon the nuclear liquid vapor phase transition. 
Because large nuclei (A > 60) are metastable objects, phases, phase coexistence, and phase tran- 
sitions cannot be defined with any generality and the analogy to liquid vapor is ill-posed for these 
heavy systems. However, it is possible to account for the Coulomb interaction in the decay rates 
and obtain the coexistence phase diagram for the corresponding uncharged system. 

PACS numbers: 21.10.Ft, 21.65.+I, 24.10.Pa 
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Since the birth of the liquid drop model, which oc- 
curred more than 60 years ago, nuclei have been recog- 
nized as charged drops of a van der Waals like fluid. Soon 
after, the concept of cold uncharged, symmetric nuclear 
matter was introduced. The characterization of its prop- 
erties, such as its phase diagram and equation of state has 
been and remains still perhaps the most eminent goal of 
nuclear physics. 

The experimental characterization of cold nuclear mat- 
ter began by setting the surface, symmetry, and Coulomb 
terms of the liquid drop expression to zero and retaining 
just the volume term. This, together with the indepen- 
dent measurement of nuclear radii (already inferable from 
the surface and Coulomb coefficients), defined the fun- 
damental properties of cold symmetric nuclear matter, 
namely its binding energy and density at saturation. 

Thus, the Coulomb interaction was reasonably elimi- 
nated from the picture in order to dispose of what was 
perceived as a troublesome inessential divergence, while, 
for better or worse, it remains all pervasive in the exper- 
imental realm. 

The experimental extension to higher temperatures 
was hampered both by the lack of a suitable container 
which became necessary to accommodate the vapor phase 
(shown to exist by the the soon discovered neutron and 
proton evaporation) and by the not obviously generaliz- 
ablc finite size effects. 

In this paper we consider the problem of the Coulomb 
interaction and the problem of the container, and we 
present a simple and natural solution to both. This so- 
lution also gives a good indication of how to deal with 
finite size with maximum generality. 

The nuclear liquid vapor phase transition, or more 
properly, the liquid vapor coexistence is being discussed 
intensely in literature QHHIllllEllli- Theoretically, 
a great deal of attention has been given to the effects of 
the nuclear finite size 0. Experimentally, discovery and 
characterization has been claimed from various quarters 
ESS US 00, although doubts remain regarding the 
validity of some approaches and the consistency of the 
results. 

Many of the theoretical approaches have been based 



upon numerical simulations of finite lattice systems. 
From these studies negative heat capacities have been 
claimed to be the signal of the phase transitions in nuclei 
and other finite systems 0. From our part we found it 
more productive to use a simple extension of thermody- 
namics which incorporates the finiteness of the nuclear 
system through the surface and other liquid drop terms 



We showed that a negative heat capacity can be triv- 
ially expected in terms of the change of the droplet (nu- 
clear) binding energy with the decreasing size of the evap- 
orating droplet (nucleus). This can be seen most clearly 
in the case of an isobaric transition: 
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Since from the Clapeyron equation 
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we can obtain the desired result for the isobaric depen- 
dence of the transition temperature on the size of the 
system 
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In all the above, A is the size of the system, AH m is 
the molar enthalpy of vaporization (AH m B(A) + T), 
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FIG. 1: The droplet dependent caloric curve at constant 
pressure. The solid line shows the drop's caloric curve 
and its dependence on radius. The scaled values of tem- 
perature, enthalpy, and radius are T' = T/AH^, H' — 
H(AH^/3c s Vi) 3 /(p_/Vi) and r' = r AH° m / (ic a VL) , re- 
spectively (see ref. pH). AH^ is the bulk molar enthalpy 
of vaporization, c s is the surface energy coefficient, and is 
the molar volume of the liquid. 



and B(A) is the binding energy per particle of a cluster 
of size A. 

In a droplet of a van der Waals liquid, AH m increases 
with increasing size A and saturates for the infinite sys- 
tems. This is of course due to the droplet molar sur- 
face energy which decreases and asymptotically vanishes 
with increasing droplet size. Thus Eq. © says that, at 
fixed pressure p the coexistence temperature decreases as 
the droplet (nucleus) evaporates. This implies a sloping 
downwards of the caloric curve in the transition region 
and a negative heat capacity in the same region. An 
example of the droplet size dependent heat capacity is 
given in Fig. ^ as the slope in the temperature-enthalpy 
correlation at constant pressure for different drop sizes. 

However, as we pointed out elsewhere js|, in nuclei 
surface effects are not the only ones to be considered. 
Coulomb effects become progressively more important 
with increasing A, Z and eventually, at about A— 60 they 
reverse the surface trend prevailing at small A values. 
Above A « 60 where the binding is at a maximum the 
binding and thus AH m progressively decreases with in- 
creasing A. The inescapable conclusion is that, within 
the scope of the Clapeyron equation, no negative heat 
capacities should be observed for A > 60, in contradic- 
tion with claims to the contrary [||. 

The above result is predicated upon the dominant 



emission of monomers (neutrons, protons), or, at any 
rate, very small clusters. Should the nucleus statistically 
breakup into large fragments, the A ps 60 upper limit 
for the observation of negative heat capacities would be 
correspondingly displaced upwards. However, not only 
are large fragments rare at low temperature but they 
never dominate at any temperature (below the critical 
temperature). In Ising models, finite or not, the average 
cluster size in the "gas phase" hardly rises above unity. 
The same holds true in experimental data 0, ■ It truly 
seems impossible to escape these thermodynamic conclu- 
sions. 

However, one may question whether the role of the 
Coulomb interaction is merely that of decreasing the 
binding energy. The long range nature of this force may 
compel us to analyze its role in more detail in first order 
phase transitions. As will be shown below, the problems 
are quite serious and threaten our ability to define a true 
first order phase transition with any generality in the 
presence of such a force. 

Let us begin with a premise. In typical first order phase 
"coexistence" , the two phases (liquid and vapor) do not 
have to be in actual physical contact. The contact sur- 
face is irrelevant as is the short range interaction between 
the two phases. In other words the equilibrium regime 
depends exclusively on the properties of each of the two 
phases, as if the other were not there. Equilibrium exists 
when the chemical potentials of the two phases are equal, 
be they in contact or not. 

Let us now introduce the Coulomb interaction in the 
problem of a drop and its vapor. 

The Coulomb interaction can be split into three parts: 
1) the drop self energy; 2) the drop-vapor interaction 
energy; and 3) the vapor self energy. 

The drop self energy, for a finite bound or metastable 
drop, is easily calculated and does not constitute a prob- 
lem. 

For the drop-vapor interaction, we consider a probe 
cluster which we can carry from the interior of the drop to 
infinity. The potential energy experienced in the process 
depends upon the particle's charge/mass and is shown 
schematically in Fig. [21 

If the particle has zero charge (top panel) , a step is ob- 
served at the droplet radius equal to the particle binding 
energy. For charges greater than zero, a maximum B c is 
observed at the approximate distance of the two droplets 
in contact. From there the potential decreases accord- 
ing to the Coulomb law and settles down at infinity to a 
value equal to the binding energy of the particle, Q. 

In this case, where we assume that any particle of any 
size is bound (Q < 0) and we forget about problem 3), 
there is no difficulty in defining a gas phase in equilib- 
rium with the droplet at infinity constituted by particles 
of all sizes whose abundance is controlled by the respec- 
tive binding energies in the standard way. The interven- 
ing Coulomb barrier B c does not alter the equilibrium, 
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FIG. 2: Top panel: the schematic potential of an uncharged 
and bound particle (Q < 0) leaving a nucleus. Bottom panel: 
the schematic potential of a bound charged particle. The 
charged particle must overcome a Coulomb barrier B c in order 
to leave the nucleus. 



although it may slow its achievement. In this case the 
vapor is constituted mainly of monomers and the coexis- 
tence pressure described by the Clapeyron equation 

dp _ Ag m , . 

dT ~ TAV m [ ' 

with the molar enthalpy AH m suitably accounting for 
both surface and Coulomb terms, is completely adequate 
to describe the liquid to vapor transition and coexistence 
(AV m is the molar volume). 

Let us now consider the case in which the probe parti- 
cle becomes unbound to the droplet above some Z value, 
due to the Coulomb interaction. Now the situation be- 
comes as depicted in Fig. In this case the droplet 
is not stable and the ground state of the system may 
consist of two or more pieces of the original drop at in- 
finity. This is naturally true already at T = 0. Thus it 
is not possible to speak properly of this drop in statis- 
tical equilibrium with its vapor, since the drop itself is 
metastable. For a nucleus like gold, the ground state is at 
least as complicated as three fragments of approximately 
size 60 nucleons at infinity. This "true" ground state is 
hundreds of MeV below the mass of the gold nucleus. 
In any statistical calculation, at any reasonable temper- 
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FIG. 3: The schematic potential of an unbound charged par- 
ticle. 



ature, one can expect a liquid-like phase consisting of a 
configuration similar to the true ground state in equilib- 
rium with some vapor. A metastable gold-like drop is an 
immensely improbable configuration because of the great 
energy chasm mentioned above. The probability of such 
a configuration can be surmised from the Boltzmann fac- 
tor P = cxp(—AE/T) where AE is the energy difference 
between the metastable state and the ground state. Es- 
timating AE ?» 135 MeV and a temperature of 4 MeV 
we obtain P ps e~ 34 or approximately 2 x 1CP 15 . 

One might argue that our point is made from energetic 
rather than free energy considerations and that it may in 
fact be incorrect. After all, equilibria work both ways, 
and typically one of the phases is at a lower energy than 
the other. 

Let us consider, then, the transition from a condensed 
phase (liquid- like) to a dilute phase (vapor-like). For an 
infinitesimal isothermal transfer, the variation of the free 
energy must be zero 

AF = AE - TAS = 0. (8) 

As we go from liquid to vapor, AE > for a typical fluid, 
but this energy increase is compensated by an equivalent 
increase in entropy, due just to the increase in molar 
volume. 

However, if AE is negative, due to the Coulomb effect, 
we need a decrease in entropy which is hardly compatible 
with expansion. 

The conclusion is that a statistical equilibrium first 
order phase coexistence and phase transition is not de- 
finable for any droplet that has unbound channels. Of 
course, the transition of the metastable droplet to its 
"true" complex ground state does not qualify as a statis- 
tical phase transition. 

This Coulomb effect seems truly devastating since it 
does not allow one to define nuclear phase transitions 
much above A sa 60. 

However there may be a solution to this difficulty. If we 
consider the emission of particles with a sizable charge, 
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we notice that a large Coulomb barrier B c is present. 
For T « B c these channels may be considered effec- 
tively closed. Consequently the unbound channels may 
not play a role on a suitably short time scale. Then a 
phase transition may still be definable in an approximate 
way. But, of course, we reach again the previous conclu- 
sion that for A > 60 heat capacities must be positive 
and therefore claims of negative heat capacities for large 
nuclei find here a most serious objection. 

Let us consider now part 3) of the Coulomb energy, 
namely the vapor self energy. As we said above, it di- 
verges for an infinite amount of vapor. For a dilute vapor, 
we could consider a small portion such that the intrinsic 
self energy /nucleon is much less than the temperature 
T. Alternatively, we could consider a finite box contain- 
ing a finite system. Unfortunately, at any other distance 
smaller than infinity the result depends annoyingly on 
the size (and shape!) of the container and on whether 
the drop is confined or not in a specified location of the 
container; a rather inelegant and non-general situation 
leading to confusing questions about true equilibrium. In 
any case, it is clear that overall, the Coulomb term makes 
the definition of phase coexistence and phase transition 
intractable and ill-posed. 

A solution to these difficulties can be arrived at only 
by asking a slightly different question: is there a way 
to obtain experimentally the signal and characterization 
of the phase diagram (transition) of a nucleus as if the 
Coulomb interaction were not there? 

As mentioned above, any attempt to define and charac- 
terize both phases in the presence of the Coulomb inter- 
action depends (at the very least) on the shape and size 
of a confining volume applied from without. This seems 
artificial and lacks a desirable generality. But nature ac- 
tually provides this "confining volume" for us. Any par- 
ticle trying to leave the nucleus is "boxed in" by a barrier 
(B s ) which depends on the particle under consideration 
and on the residual nucleus (or the "complement" ) . The 
top of this potential barrier is close in shape to the po- 
tential of two objects, particle and complement, in near 
contact. The tops of these barriers are actually condi- 
tional saddle points 0], conditional in the sense that 
the mass asymmetry is considered frozen. 

According to standard transition state theory all these 
saddles are in statistical equilibrium with the droplet and 
the decay rates give direct information on their popula- 
tion which is naturally controlled by a Boltzmann factor 
exp(— B s /T). In particular for large enough B s the ob- 
served experimental abundances are directly related to 
first chance emission and thus to the transition state rate. 

Now B B is composed by 



B s — Et 



rfacc 



E 9 ' 



■face 



E, 



Coul 



(9) 



where -B s s urface and E^ ifacc are the surface energies of 
the saddle and ground state respectively, and E^ oul and 
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FIG. 4: A schematic representation of the Coulomb correc- 
tion when the emitted fragment is bound (left panels) and 
unbound (right panels). In both cases one can remove the 
Coulomb energy of the saddle configuration and calculate the 
Q value using surface energies only (bottom panels). The re- 
sulting hypothetical gas will be composed of fragments that 
are bound to the droplet (Q sur f aC c < 0) for all fragment par- 
titions. 
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FIG. 5: The dependence of the saddle barrier (B s , dotted line) 
is plotted as a function atomic number for a schematic binding 
energy which includes only Coulomb and surface terms. Each 
saddle barrier consist of two parts: the — Q value (short dash) 
and the Coulomb barrier B c (long dash). After correction for 
Coulomb, the relevant Q value calculated using surface terms 
only binds the fragment to the drop and is shown by the solid 
line. This is calculated for a Z=50, ^4=118 nucleus. 



-^Coui are tne Coulomb energies for the same two con- 
figurations. Since the Coulomb energies can be easily 
estimated assuming a two touching spheres configura- 
tion for the saddle and one sphere configuration for the 
droplet (see Fig. 0J, we can correct the rates by divid- 
ing away the Boltzmann factor containing the Coulomb 
terms and be left with only the rates/ abundances per- 
taining to the decay of an uncharged drop, for which all 
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channels are bound by the extra surface energy Qsurface = 
Surface ~ Surface ( see Fig.0. These rates are now inde- 
pendent of distance and are proportional to the effective 
partial concentration of the hypothetical gas in equilib- 
rium. We speak of a virtual gas phase because it is not 
and it need not be present. This picture of a free evap- 
oration of a droplet in vacuum neatly bypasses also the 
need for a physical presence of the vapor. The result- 
ing situation is very much that described by the Fisher 
droplet model for the composition of a saturated va- 
por in equilibrium with a liquid droplet. The Fisher 
droplet model can be directly co-opted to describe the 
(first chance) fragment abundances of a nuclear physics 
experiment after correction for Coulomb effects. From it, 
it is easy to obtain the coexistence diagram for any nu- 
clear system deprived of the Coulomb interaction |]J, |2j . 
This is in the same spirit as in nuclear matter calcula- 
tions in which neutrons and protons are considered as 
distinct particles, but without any Coulomb interaction. 

In summary, we have illustrated the difficulties brought 
about by the Coulomb interaction in the problem of the 
nuclear heat capacities. We have shown that for systems 
with bound decay channels, negative specific heats, if 
definable at all, are allowed only for A < 60. For unbound 
channels, namely nuclei metastable to breakup into two 
or more large fragments, no phase description or phase 
coexistence with the normal thermodynamic generality is 
possible. However, we have shown that decay transition 
rates can be corrected for the Coulomb interaction in a 
natural way and the phase diagram for the corresponding 
chargeless system can be easily obtained. 
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